Superconductivity assisted by inter-layer pair hopping in multi-layered cuprates 
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In order to explore why the multi-layered cuprates have such high 7* c 's, we have examined various inter-layer 
processes. Since the inter-layer one-electron hopping has little effects on the band structure, we turn to the inter- 
layer pair hopping. The superconductivity in a double-layer Hubbard model with and without the inter-layer 
pair hopping, as studied by solving the Eliashberg equation with the fluctuation exchange approximation, reveals 
that the inter-layer pair hopping acts to increase the pairing interaction and the self-energy simultaneously, but 
that the former effect supersedes the latter and enhances the superconductivity. The inter-layer pair hopping 
considered here is for off-site pairs, for which we discuss the effect of retaining SU(2) symmetry, along with 
how the sign of the pair hopping determines the relative configuration of d-waves between the adjacent layers. 

PACS numbers: 74.20.-z, 74.62.-c, 74.72.-h 



Introduction. — Although we are witnessing the discovery 
of new classes of superconductors that include the iron-based 
and organic superconductors|uJ|2|], the high-Tc cuprate super- 
conductors stand out in having the highest- T c to date. Specif- 
ically, among various families of the cuprate, the highest T c 
occurs in the multi-layered cuprates that have CUO2 planes 
in a unit cell, typically the Hg-series HgBa2Ca n _iCu„02 n +2+a> 
where T c depends on the number, n, of the CuC>2 planes with 
T c increasing for n = 1 to 3 and decreasing slightly for n > 4. 
is the highest T c superconductor[3|. Experimentally, the elec- 
tronic band structure has been probed with ARPES for the 
Bi-based triple-layered cuprate (Bi-2223)04j. Another ex- 
periment examines the optical Josephson plasma modes aris- 
ing from inter-layer Josephson couplings from the reflectivity 
spectra in the Hg-based multi-layered cuprates for n = 2 - 5, 
where the change in the Josephson coupling strength is shown 
to be correlated with T c |5[] . 

There have been several theoretical studies for multi- 
layered cuprates: Refs. [a, lTJ discuss that an inter-layer pah- 
hopping that arises as a process second-order in inter-layer 
one-electron hopping enhances the superconductivity, Ref. 
[8] treats the inter-layer pair hopping macroscopically in a 
Ginzburg-Landau scheme, Ref. H examines Coulomb en- 
ergy in layered structures, and Ref. [10] looks at an effect of 
inter-layer one-electron hopping for a double-layer Hubbard 
and t-J models. Given the situation, our purpose here is to mi- 
croscopically investigate mechanism of the superconductivity 
in multi-layered cuprates focusing on the effects of inter-layer 
pair hopping. 

Motivated by this, here we start from a double-layer 
Hubbard model to explore microscopically the multi-layered 
cuprates by examining various inter-layer processes. The 
inter-layer one-electron hopping has turned out to exert lit- 
tle effects on the first-principles band structure (not shown), 
so that we turn to the inter-layer pair hopping. The hopping of 
Cooper pairs across the layers should in general exist as a ma- 
trix element of the long-range Coulomb interaction|)TT 121. 



and this should affect the superconductivity as a process in- 
trinsic in multi-layer systems, but whether and how the super- 
conductivity is enhanced has not been well understood. Since 
we are talking about d-wave pairing that basically mediated 
by antiferromagnetic spin fluctuations around specific regions 
in fc-space, we have to adopt a method that can incorporate k- 
dependent pairing interactions. Hence we adopt here the fluc- 
tuation exchange approximation (FLEX) lll3l4l7ll . whose re- 
sult is fed into the Eliashberg equation. We shall show that the 
inter-layer pair hopping acts both ways to increase the pairing 
interaction and decrease the quasi-particle life time (with an 
increased self-energy), but the former effect is found to su- 
persedes the latter and enhances the superconductivity. The 
inter-layer pair hopping considered here is for off-site pairs as 
necessitated if we want to treat d-waves, for which we discuss 
the effect of retaining SU(2) symmetry, along with how the 
sign of the pair hopping determines the relative configuration 
of d-wave between the adjacent layers. 

Formalism. — We consider a doule-layer Hamiltonian H 
with the inter-layer pair hopping ff pa ir. 



H = H, + H v + H, 



1 pair » 



(l) 



where the one-electron kinetic energy, H t = 
~2Z a /3 2Zij}Za- t "f c 'i^ c ^ a -' an d the Hubbard interaction, 
Hu — UYja 2ji c ? c if c< ii c %' ^ defined in a usual way with 
creating an electron at i-th site with spin cr in the layer a, 
t.j the transfer integral and U the on-site Coulomb repulsion. 
H t consists of the intra-layer (a = p) and inter-layer (a + /3) 
one-electron hoppings, where the intra-layer component is 
here considered for the nearest-neighbor t = 0.5eV, second- 
neighbor f = O.leV up to the third-neighbor t" = 0.08eV. 
For the inter-layer one-electron hopping we take a usually 

adopted forrrL ft, = 2<^g Zk Za- 1 (cos k x - cos k y ) cQc/j^ 
in fc-space lll8[ Il9ll . with t z = 0.05eV here. These values 
of the one-electron hoppings are basically determined by 
a downfolding from the first-principles bands, but here we 
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make a simplification in which we take common values 
between the single- and double-layer cases for a transparent 
comparison. 

Now the question is the form of the inter-layer pair hopping 
//p a i r . Here we take a rather general form, 



off 

pair ±± pair ' ^pair' 



#naii - H°*- + H' 



(2) 



where, on top of the usually considered inter-layer on-site pair 
hopping, 



"parr U Zj Zj <T 4 4 fl" 

ct+p i 

we also consider inter-layer off-site pair hopping, 
H° B = H° m) + H° B(2) 

pair pair pair ' 



(3) 



(4) 



where the first term, 

nn 

pair iT jl jl if v ' 

a±P ij 

is the hopping of a spin-singlet pair formed on nearest- 
neighbor intra-layer sites from one layer to another, with 
denoting a sum over nearest-neighbors. In addition, we have 
to note that, if we want to preserve the spin SU(2) symmetry, 
we should include 



H 



off(2) 
'pair 



U' 



(6) 



in which the spins of the pair are exchanged during the hop 
(FIG.Q]). While the on-site term is considered to be the largest 
inter-layer pair hopping, the hopping of off-site pairs should 
be not only the second largest inter-layer pair hopping aris- 
ing from long-range Coulomb interaction, but may also play 
a crucial role, since we are talking about an off-site, d-wave 
pairing. 




FIG. 1: (color online), (a) Schematic inter-layer hopping of on-site 
pairs (if™ ) and off-site pairs (# off ). The latter consists of H° m> 

r v pair ' t v pair' pair 

and H °? 2 \ where the spins of the pair are exchanged during the hop 
in H°^\ (b) Diagrams for the non-spin-flip interaction H^r (left 
panel) and spin-flip interaction H°P ' (right). 



Now, the FLEX approximation, which is a conserved ap- 
proximation with bubble and laddar diagrams included lll3l 



11411 . is one of the standard methods for treating the spin- and 
charge-fluctuation mediated pairng self-consistently with the 



self-energy effect incorpratedlfl5 -171. Let us start with show- 
ing that the method can be extended for treating the pair- 
hopping processes introduced here. Derived from Dyson- 
Gor'kov equation, the linearized Eliashberg equation for the 
gap function A a p(k) reads, in the present case, 



H k' a'P' y5 

x G a , y (k')A y6 {k')G p > s {-k'). 



(7) 



Here k = (k, cl>„) is the two-dimensional wave number and 
Matsubara frequency for fermions, f3 = l/T {kg = 1) and A the 
eigenvalue of the Eliashberg equation, where Tq is identified 
from A — 1 but A also serves as a measure of the strength of 
superconductivity. The pairing interaction in the present case, 
being involved with layer index, becomes a bit complicated (a 
2x2x2x2 tensor) as 



V pan " (a) 



u + l- 



UXoU 1 UxoU 



(q), (8) 



where \xo] aa >0p(q) = -(l/N/3)ZkG(i a (k + q)G a ,p(k) is the 
polarization function, while U, also a 2 4 tensor, represents 
the interaction, which can be expressed as a 4 x 4 matrix, 
with the four rows (columns) corresponding to aa'(J3/3') = 
11,22, 12, 21, as 



U(q) 



(U 














u 

















U' + U"(q) 


lo 





V + U"(q) 


J 



(9) 



with U"(q) = 2U" (cos q x + cos ^ v ). 

For each layer, the d-wave pairing is favored by the intra- 
layer pairing interaction ( V^aaaiq) m the p resent notation) that 
has peaks around Q - (+7r. +7r) lll5l4l7l 12011 . For a double- 
layer model having inter-layer pair hoppings the question is 
how the inter-layer pairing interaction V^ a (q) (a + f3) affects 
superconductivity. 

Results. — Now we present the results comparing the situ- 
ations in the presence and absence of the inter-layer pair hop- 
ping in FIG. 12 This plots the eigenvalues of the Eliashberg 
equation A against the band filling «, where we set U = 2.5eV 
here, a relatively small value compared to the realistic param- 
eter but appropriate to FLEX which is a weak-coupling for- 
malism. For the pair-hopping interactions, we set U' = 0.5eV 
and 2U" = -0.5eV, values chosen to be much smaller than 
U but still significant, while the effect of the sign U" will be 
discussed later. 

Beside the Eliashberg A, we also display in FIG.[3]the inter- 
layer pairing interaction Vf^(q) (= V^ 2 (q)) (at n = 0.85). 
This is important, since the d-wave pairing within each layer 
has a strongly ^-dependent form, An(£) = A22(k) ~ cos k x - 
cos k y , so that the real question, for multi-layer cases, should 
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FIG. 2: (color online). The eigenvalue A of Eliashberg equation 
against the band filling n for the double-layer model. Black (dashed) 
line: no inter-layer pair hopping, purple: with //™ r only, green: with 

and 



H° mi} only, blue: with H oa + If 

pair ' pair pair 

isolated diagrams for // off<1) 



off! 1 ) 



-off(l) 



, pink: with H on + H' 

' r pan- pair 



of// om "andtf' 

pair pair 

mesh for (k x , k y , a>„) 



, red: with all of _/¥ p!i ;, except for mixing 
O.OleV with a 32 x 32 x 2048 



pall 

Here we set T 




FIG. 3: (color online). Real part of Vf 221 (q, £„=o)) when various pair- 



ing interactions are switched on one by one at n = 0.85. (a) ff™ ir 

, (d) H°". + H off{>) and isolated 

pair J ' pair pair v ' pair pair 

diagrams for H°^\ (e) All of // pa j, except for mixing of //° ff<1 1 



only, (b) H° m) only, (c) H on + H° m) 

J 7 v ' Dair J ' v ' pair p; 

;) , (e) All of H„ 



H' 



i>ll[2) 



pair 



and 



be the effect of inter-layer pair hoppings on such anisotropic 
gap functions. 

In order to resolve the effects from various terms, let us 
switch on the terms one by one. First, black dashed line in 
FIG. |2]represents the result of the double-layer model without 
inter-layer pair hopping. FLEX becomes unreliable when the 
band filling becomes too close to the half-filling, so that we 
only plot the result up to n < 0.9. Now, purple line in FIG. 



|2] is for the model with inter-layer on-site pair hopping H™ r 
only, where the pair hopping is seen to suppress the supercon- 
ductivity in fact. This result, which may at first seem strange 
since an inter-layer pairing interaction would naively enhance 
superconductivity, comes from the following fact. The inter- 
layer pair hopping does produce an inter-layer pairing interac- 
tion as displayed in Fig. 3(a), which is expected to enhance the 
intra-layer superconducting gap functions An(£) and A22W 
in the sense of Suhl-Kondo mechanism. However, the inter- 
layer pair hopping also increases the (intra-layer) self-energy. 
An increased self-energy is a bad news for superconductiv- 
ity, and the result here indicates that this effect supersedes the 
enhanced inter-layer pairing interaction. If we look at FIG. 
Oa), the inter-layer pairing interaction V^^iq) only shows 
barely visible peaks around Q. This is because the on-site 
inter-layer pair hopping Hamiltonian has no ^-dependence to 
start with, and FLEX diagrams do not render a significant k- 
dependence. This is why V\ Z7A (q) is insufficient for overcom- 
ing the increased self-energy. 

By sharp contrast, green line in FIG. |2] which represents 
the result when one of inter-layer off-site pair hopping, H°^\ 
is switched on (without //™ ir ), exhibits a significant enhance- 
ment. Indeed, FIG. 0b) shows that V^ 22l (q) develops a signif- 
icant ^-dependence, and this is how the enhanced pairing in- 
teraction overcomes the increased self-energy, since H°^\cc 
cosq x + cosq y ) originally possesses a large ^-dependence, 
where the peaks at Q are intensified by the pairing interaction. 

Now the question is whether the addition of the on-site pair 
hopping (fKj) degrades the enhancement due to H°^\ Blue 
line in FIG. |2] representing this situation shows that the su- 
perconductivity is enhanced even above the case when 
alone is switched on. This may first seem to contradict with 
the fact that ff™ r suppresses the superconductivity, but, if 

off(l) 
pair "pair 
off 



we go to FIG. He), V^(q) with both of and H' 

switched on is more reinforced around Q than when H' 

1 

alone is present. The increase in V^^iq) is caused by the 
process in which //™ ir raises the peaks of V^j (q) around Q 
from // off(1) through the spin-fluctuation term I Ux £ u in ©. 

pair 2 1-Uxo 

Let us finally discuss the effects of the other inter-layer 



off-site pair hopping H 



-off(2) 



The term is required for the 



SU(2), but, being a spin-flip interaction, does complicate 
the diagrams as follows. The term reads in £-space as 
<r 2> = (-VN)Z^Z kk , q U"(pc? +ql c$ + as de- 

picted in FIG. |TJb). We can readily extend the FLEX when 
all the interactions are of the spin-flip form, where the formu- 
lation is similar to the usual FLEX. 
Therefore we first take account of H° 



+ ^off( l) an( j j sQ _ 



lated diagrams for separately (i.e., excluding the mix- 



ing of H' 



off(2) 
pair 



with spin-nonflip H on + // off(1) ). In this case, 

r r pair pair 7 



the pairing interaction © is replaced with W du [U, U', U"] + 
W dir [0, 0, U"]. For details, see the supplemantal material lT2lll. 
Dramatically, the addition of // off<2) is seen as pink line in 

J pair r 

FIG. |2] to enhance the superconductivity much more than the 
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case with H° n and // off(I) alone. A reinforced V v ,fL (q) around 

pair pair 1 22 1 v ^ ' 

Q are indeed seen in FIG. Ed). 

Finaly, we take account of the mixing of //™ r + H^Jf with 

^° ^ eat m is rigirously is difficult because they have 
respective ^-dependences, and their mixing acts as a kind of 
vertex corrections (see FIG. 1 of supplemental material[21 j). 
However, we have confirmed from the selfenergy that the ef- 
fect of the vertex corrections is numerically negligible, so 
that we can take account of all of // pa ir except for the mixing 
of H°2 n and by replacing © with V p ' dir [U, U', U"} + 

yP air [f/, U', U"] - yP air [t/, U', 0]. The first (second) terms rep- 
resent the non-spin-flip (spin-flip) interactions, while the third 
term subtracts the double counting[21]. Red line in FIG. |2] 
represents the result in this scheme, where the superconduc- 
tivity is enhanced even above the pink one. FIG. [3je) con- 
firms that V^iiq) i s more reinfored around Q than in FIG. 
[3jd). The increase in V^^iq) (pink line to red in FIG. |2} is 
caused in FLEX because a combined effect of H on and H° s ^ 2) 

pair pair 

raises V^t (q), as a combined effect of H° n and H°. raises 

1221 pair pair 

V^22i(l) (green line to blue). 

Now, we are in position to construct a phase diagram of 
the double-layer system, in which we can compare the re- 
sult with and without inter-layer pair hopping H pa j, in FIG. 
|U Superconducting (SC) phase boundary is identified from 
the eigenvalue of linearized Eliashberg equation A reaching 
unity. The antiferromagnetic (AF) phase boundary is deter- 
mined in a usually adopted way from the (in the present case 
the intra-layer) [UxfAaaaa approaching unity (0.975 here). 



AT C ~ 0.02 x t ~ IQQK. On the other hand, AF transition tem- 
perature Tn for the double-layer model with inter-layer pah- 
hopping // pa j r slightly decreases from the case without, which 
is because the divergence of the spin susceptibility Xaaaa * s 
suppressed by the self-energy increased due to the inter-layer 
pair hopping. 

Let us finally discuss the effect of the sign of U" in H°^ k . 
We can show, as seen in FIG. |5(a)] that the sign is important 
in determining the relative configuration of the An and A22- 
The inter-layer pairing interaction Vj^j (q) with U" < 0, as 
we have assumed so far, favors the configuration where the 
anisotropic, in-plane gap functions, An and A22, are arrayed 
in-phase, which is dictated by inter-layer pairing interaction 
in FIG. [3] If we change the sign to U" > 0, we end up with 
a configuration where An and A22 are arrayed out-of -phase 
(i.e., rotated by 90 degrees with each other) as shown in FIG. 
|5(b)| In this case, Vf^iq) also changes sign. 




(a)C/" < 0. (b)U" > 0. 



FIG. 5: (color online). Configuration of the d-wave, in-plane gap 
functions, An, A 2 2 for different signs of U" 




1.00 



FIG. 4: (color online). Phase diagram on T and n (carrier concentra- 
tion) for the double-layer system with (red lines) and without (black) 
inter-layer pair hopping H v ^ n . T c is SC transition temperature while 
7" N AF transition (Neel) temperature. The arrow represents the in- 
crease of T c arising from the inter-layer pair hopping. 



As is seen in FIG. |4j SC transition temperature T c for 
double-layer model in the presence of the inter-layer pair 
hopping //p a i r is higher than the case in the absence for all 
the range of the carrier concentration considered here. For 
U' = 0.5eV and 2U" = -0.5eV, the increase of T c amounts to 



To be more precise, however, the configuration is not de- 
termined solely by the sign of U": even in the absence of the 
inter-layer pair hopping, the in-phase configuration is favored 
through the off-diaginal Green's functions, Gyi and G%\, in 
the Eliashberg equation (0. When the inter-layer pair hop- 
ping is switched on, the effect of V^ 22l (q) has to overcome 
this effect of V^j(^) which favors the in-phase configuration 
before out-of-phase configuration is realized for large enough 
U" > 0. 

Summary. — To summarize, superconductivity in a double- 
layer Hubbard model with and without the inter-layer pair 
hopping is studied by solving the Eliashberg equation with 
the fluctuation exchange approximation. We have shown that 
the inter-layer pair hopping acts to increase both the pairing 
interaction and the self-energy, but that the former effect su- 
persedes the latter and enhances the superconductivity. The 
inter-layer pair hopping considered here is for off-site pairs, 
for which we have found that the extra off-site pair-hopping 
term needed to preserve SU(2) symmetry actually acts to en- 
hance the superconductivity even further. We then end up 
with a phase diagram for the double-layer model where the su- 
perconducting boundary is significantly higher than the case 
without inter-layer pair hopping. 

In evaluating the present mechanism, an estimate (e.g., with 
c-RPA) of the magnitude of inter-layer off-site pair hopping 



ff°2j in rea l materials should be important. One possibly rel- 
evant quantity is the optical Josephson plasma modes, which 
have been observed in ref [5] for Hg-based cuprates with 2- 
5 layers, where the mode energy of the double-layer system 
is ~ 800 cirT 1 that corresponds to ~ O.leV. It is an inter- 
esting future problem to examine the actual relation of this to 
the inter-layer pair hopping considered here. Also, the case of 
three or larger number of layers is interesting, for which the 
study is under way. 
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PRESTO. 



[1] H. Aoki, J. Supercond. Novel Magnetism 25, 1243 (2012). 
[2] Y. J. Uemura, Nature Materials 8, 253 (2009). 
[3] A. Schilling, M. Cantoni, J. D. Guo, and H. R. Ott, Nature 363, 
56 (1993). 

[4] S. Ideta, K. Takashima, M. Hashimoto, T. Yoshida, A. Fuji- 
mori, H. Anzai, T. Fujita, Y. Nakashima, A. Ino, M. Arita, 
H. Namatame, M. Taniguchi, K. Ono, M. Kubota, D. H. Lu, 



Z.-X. Shen, K. M. Kojima, and S. Uchida, Phys. Rev. Lett. 
104, 227001 (2010). 

[5] Y. Hirata, K. M. Kojima, M. Ishikado, S. Uchida, A. Iyo, 
H. Eisaki, and S. Tajima, Phys. Rev. B 85, 054501 (2012). 

[6] P. W. Anderson, Theory of Superconductivity in the High-Tc 
Cuprate Superconductors (Princeton University Press, 1997). 

[7] S. Chakravarty, A. Sud0, P. W. Anderson, and S. Strong, Sci- 
ence 261, 337 (1993). 

[8] S. Chakravarty, H. Kee, and Volker., Nature 428, 53 (2004). 

[9] A. J. Leggett, Phys. Rev. Lett. 83, 392 (1999). 
[10] S. Okamoto and T. A. Maier, Phys. Rev. Lett. 101, 156401 
(2008). 

[11] K. Kusakabe, J. Phys. Chem. Solid 73, 1546 (2012). 
[12] K. Kusakabe, J. Phys. Soc. Jpn. 78, 114716 (2009). 
[13] G. Baym and L. P. Kadanoff, Phys. Rev. 124, 287 (1961). 
[14] G. Baym, Phys. Rev. 127, 1391 (1962). 

[15] N. E. Bickers, D. J. Scalapino, and S. R. White, Phys. Rev. 
Lett. 62, 961 (1989). 

[16] N. Bickers and D. Scalapino, Ann. Phys 193, 206 (1989). 

[17] T. Dahm and L. Tewordt, Phys. Rev. Lett. 74, 793 (1995). 

[18] O. K. Andersen, O. Jepsen, A. I. Liechtenstein, and 1. 1. Mazin, 
Phys. Rev. B 49, 4145 (1994). 

[19] O. K. Andersen, A. I. Liechtenstein, O. Jepsen, and F. Paulsen, 
J. Phys. Chem. Solid 56, 1573 (1995). 

[20] K. Kuroki and H. Aoki, Phys. Rev. B 56, R14287 (1997). 

[21] See Supplemental Material for the FLEX with intra- and inter- 
layer interactions in a two-layer model. 



1 



Supplementary information for 
'Superconductivity assisted by inter-layer pair hopping in multi-layered cuprates" 



We present an outline of the extension of the FLEX (fluctu- 
ation exchange approximation) to include the spin-flip as well 
as non-spin-flip interactions. 

Multi-orbital FLEX — In the present context this we start 
with re-formulating the FLEX with intra- and inter-layer inter- 
actions in a two-layer model. We first separate non-spin-flip 
interactions such as H v , H° 



actions such as // off(2) . The non-spin-flip part (Hi 



and //° ff(I) from spin-flip inter- 
pair pair * r 



pair 

H°^) can be expressed as 



H° n and 

pair 



#nsf - 



Z Z u ^w€A-Atfv c 10 ) 



kk' q afifi' a' 



o 




FIG. 6: Second-order cross term between H 



off(I) 



and//' 



011(2) 



where a, fl, etc denote the layer, q the momentum trans- 
fer, and the nonzero components in the present model are 

tfffnfo) = VgnW = U, U^(q) = U£ n {q) = 
U' + U"(q). On the other hand, the spin-flip term 

tfplf = -(VN)^, q ^U%q)ct q ftAAv Which 
is required for SU(2) be preserved, can be expressed as 



lation operators, as 



it _ _ U V V of »t a a 
U ~ N Z^ C t+ql C k'-qfk'l C kV 



N 



kk' q a 



rron 
pair ft 



y y c oi c° ] ' </> a 



kk'q atp 



(12) 



= -^Z Z u UMl 

kk' q afifi'a' 



(ID 



where the form cjcjlcjCf signifies the spin-flip, and the 



nonzero components in the present model are Uf 22l (q) 



U"(q). While Hjj and //™ r can also be expressed 

off(l) 
pair 



in a spin-flip form (see below), we cannot cast both of H[ 

and simultaneously into a single expression like above 

if we want to have the prefactor as a function of q. 

In FLEX, all of the bubble and laddar diagrams composed 

of H n sf and H s f have to be summed, which include cross terms 

of // off(1) and H° S(2 \ It is difficult to treat the cross terms ex- 
pair pair 

actly, since a kind of "vertex correction" as shown in FIG. 
[6] exists already in the second-order in the perturbation ex- 
pansion. Fortunately, however, we have confirmed here nu- 
merically that such diagrams are much smaller than the other 
terms in the same order, which is partly because the momen- 
tum dependence is different between // off(1) and // off(2) . We 

r pair pair 

can therefore ignore the diagrams composed of the mixing of 
// off(1) and// off(2) . 

pair pair 

With this we can actually sum all the bubble and laddar 
diagrams for both of ff ns f and H s f, which is performed as 
follows. The FLEX for H v , H on and H° m> (i.e., all the 

u ' pair pair v ' 

bubble and laddar diagrams composed of H ns f) can be per- 
formed in a standard way, where the only difference is to take 
into account the tensorial interactions and susceptibilities (i.e., 

u(q),Mq))- 

We next take account of the mixing of Hn, H on . . with H°J^ 

" pa.ii pan 

employing the following technique. First, we cast Hjj and 
H°" h into a spin-flip form, by rearranging creation and annihi- 



7/° ff(2) in the spin-flip form have nonzero 

pair 1 1 

components Uf m (q) = Uf 222 (q) = U, Uf 22l (q) = Uf u2 (q) = 
U' + U"(q). We can now take account of the mixing of Hjj, 



The H v and H on + 

u pair pair 



^pair anc ^ ^pair 2) wnen a U of the bubble and laddar diagrams 
composed of H s f are summed (see FIG. [7). When we use 
the technique above, the effective interaction for the normal 
self-energy composed of ff nsf and H sf , V Gasf [U, U', U"(q)] 
and V Gsf [U, U', U"(q)] respectively, are equivalent, and the 
pairing interaction for the anomalous self-energy composed 
of ff nsf and H sf , V Fnsf [U, U', U"(q)] and V F sf [U, U', U"(q)] 
respectively, are also equivalent. 



Host 



cross term 




FIG. 7: FLEX can be performed for each of the Hamiltonian compo- 
nents encircled by ovals. In addition, cross terms exist between the 
components indicated by an arrow. 



Finally, the diagrams composed of H ns f and those composed 
of H s f are added, but we have of course to subtract the double- 
counted diagrams composed of Hu and H°" h .. This is achieved 

by putting the effective interaction V eff for normal self-energy 
IF and the pairing interaction V pdn for the anomalous self- 
energy S F as 



y err = ycnstr^ [/ ' ) u"(q)] + V Ust [U, U' , U"(q)] 
- V G(n)sf [t/,t/',0], 



(13) 



2 



where 



V pair = V Fnsf [f/, U', U"(q)] + V Fsf [U, U', U"(q)] 
-V F(n)si [U,U',0]. 



(14) 



We first write down the multi-orbital FLEX with intra- and 
inter-layer interactions belonging to H ns {. The normal self- 
energy for the inter-layer interactions is given as 

^ nsf » = ^ Z Z v?^(* - k')G a ,,,(n (is) 



fc' a'/ 



where 



VG.nsf 



(q) = V-^te) + V OL (q), 
nsf 



yCoB f , _ 



UxoU 



. 1 - UxoUxo 



(«), 



UxoUxqU 



l-Uxo 



nsf 



(16) 
(17) 

(18) 



a'aP'P 



Here y G oB is the bubble-diagram contribution to the effective 
interaction for the normal self-energy, where odd numbers of 
bubbles are included due to the spin selection rule in H ns f, 
while V CL is the ladder-diagram contribution to the effective 
interaction. The polarization function is defined as 



\X0]aa'/3fi'(q) 



^YjGpJk + q)G ar (k), (19) 



N/3 



V Fns{ (q) = V FeB (q) + V Fh (q), 



u 



nsf 



1 - UxoUXola'aflj 
nsf 

(«), 



(q), 



UxoU 



l-Uxo 



(24) 
(25) 

(26) 



with the same rule for the tensor products for V FeB (q) and 
V FL (q) as in the normal self-energy above. Therefore V F (q) is 
written as 



(j+ 3_UxaU 1 UxoU 



2l-Uxo 21 + Uxo 



nsf 



(?)■ (27) 



Now we turn to the multi-orbital FLEX with intra- and 
inter-layer interactions belonging to the spin-flip H sf . The nor- 
mal self-energy for the inter-layer interactions is given as 

Z C a f(k) = ± Yj Z V «'«y k - WMk')- (28) 



k' a'l 



For the spin-flip H s{ we have to take account of all of bubble 
diagrams and odd numbers of ladders due to the spin selection 
rule in H s f . 

However, we end up with the same form for the effective in- 
teraction for the self-energy V G sf (q) as before, with separated 
spin and charge fluctuation parts, as 



which isa2x2x2x2 tensor and can also be expressed as a 
4x4 matrix. As for the products of tensors, we have 



kk' 

for V G oB (q), and 



(20) 



(21) 



for V G h {q). 

For the non-spin-flip part with the on-site Hubbard inter- 
action V aaaa and the inter-layer Cooper pair hopping terms 
V a pp a (a + /J), the tensor products above are equivalent, and 
we arrive at 



| 3 UxoU | 1 UxoU 



2l-UJe 2 1 + c/^o 



UxoU 



nsf 



(q), 



ia'ap'p 



(22) 

where the second (third) term on the right-hand side is the 
spin- (charge-) fluctuation part. 

The anomalous self-energy for the inter-layer interactions 
is given as 



^; sf » = 4 Z Z K^Ak - V)F„r&), 



(23) 



f 3 UxoU , 1 UxoU fj , f7 

U + — + UxqU 

2l-Uxo 2 1 + U X o 



(?). 



a'aP'P 



(29) 

Similarly, the anomalous self-energy for the inter-layer inter- 
actions is given as 



(30) 



k' a'P' 



where we have to take account of all of bubble diagrams and 
the even number of ladders due to the spin selection rule for 

#sf- 

Thus we again end up with the same form for the pairing 
interaction for the anomalous self-energy as 



yF.sf 



(«) = 



0+ 3_UxoU 1 UxoU 



2\-U X o 2\ + Uxo\ a , m 



(?), (31) 



app' 



with the spin- and charge-fluctuation parts. 

Finally, the normal and anomalous self-energies are written 

as 



^ = 4 Z Z y «W fe - 



N/3 



k' a'p 



(32) 



k' a'l 



where V eS and V pair are expressed as Eq. (fT3l and (fl4] >. re- where A(£) = f. F (k). 
spectively. If we plug these into Dyson's equations for the 
anomalous Green's functions, we have the Eliashberg equa- 
tion, 



^ fc' »'/?' y<5 (33) 

x G^y{k'^ 5 (k')G^ s {-k'), 



